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^ ■ Abstract 



o 



We analyze the spectra of generalized Fibonacci and Fibonacci-like operators in 
Banach space I . Some of the results have application in population dynamics. 



in 

1 Introduction and preliminaries 



dcf 



Let I 1 denote the Banach space of all real sequences x = (x\,X2,x^, • • • ) such that \\x\\i 

\ x k\ < °°- Let H : I 1 — > I 1 be a linear operator on I . The resolvent set of H, p(H) is 
the set of all complex numbers A such that the operator XI — H has a bounded inverse, 
where I : I 1 — > I 1 is the identity operator. The set cr(H) = C \ p(H) is the spectrum 
of H. The spectrum is further subdivided into three mutually disjoint parts, the point 
spectrum a p (H), the continuous spectrum a c (H) and the residual spectrum a r (H). The 
point spectrum is the set of all A € C such that XI — H has no inverse. As in the finite 
dimensional case, such A are also called eigenvalues and the corresponding non-zero vectors 
x G I 1 , such that (XI — H)x = are called eigenvectors. The continuous spectrum is the 
set of all A not in p(H) or (T P (H) for which the range of XI — H is dense in I . The residual 
spectrum is the set of all A in o(H) which are not in (J P (H) or a c (H). The spectral radius 
of H is 

r CT (tf) d = f sup |A|. (1) 

Xea(H) 

The operator H has a matrix representation H in the standard basis == 6^ , where 
is the Kronecker symbol. 

We shall also use two standard results: first, if the operator H is bounded or closed and 
has a matrix representation H, then the transpose matrix H* is the matrix representation 
of the operator H* : l°° -> l°° and (see e.g. [TaLa80j . |Gol661 Corollary II.5.3] or [Hal55l 
Theorems 3.2 and 3.3]) 

o-„(ff*) Q °p(H) U a r (H), a r (H) C a p {H l ). (2) 
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Second, if H is bounded, then (see for example [TaLa8CH (3-5)]) 



r a (H) = lim \\H k \\l /k . (3) 

k— >oo 

Our aim is to classify spectra of two classes of generalized Fibonacci and Fibonacci- 
like operators. For the first class of operators their spectral radii are expressed in terms 
of largest real positive roots of certain polynomials and the coefficients of their powers 
behave like generalized Fibonacci sequences, as we shall see in section 2. 

The second class of operators, which also has applications in mathematical biology, is 
analyzed in a similar manner in section 3. 



2 Generalized Fibonacci operators 

Let the linear operator F n : I 1 — >■ I 1 be defined by 



(xi,x 2 ,x 3 ,---) ->• I 2^ x k,X!,x 2 ,x 3 ,- ■ ■ \ , n = 1,2,3,. 

\k=n+l / 

Each F n is bounded and its matrix representation in the standard basis is 



(4) 



/ ••• 1 1 111 

10 00 000 

01 00 000 

00 10 000 

00 01 000 



V 



\ 



(5) 



/ 



Following the analysis of the spectrum of F\ by Halberg |Hal55j , the spectrum of F n 
is classified in several steps which are summarized as follows: 



1. first, by solving the equation 

(XI -F n )x = 0, x^0, 
we show that the point spectrum is 

<J P (F n ) = {A e C : A n+1 - \ n - 1 = 0, |A| > 1}, 

2. second, by solving the equation 

{XI -F n )x = y, x £ 0, 



(6) 
(7) 

(8) 



we compute the inverse (XI — F n ) 1 and show that the resolvent set consists of all 
A such that |A| > 1 which are not in a p (F n ), that is, 



p(F n ) = {X g C : |A| > 1, A" +1 - X n - 1 + 0} 



(9) 
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3. third, we analyze the transposed operator and show that 



a p (Fl) = {A € C : |A| < 1, A + 1}, 
which, together with ([2]), implies that the residual spectrum of F n is 

a r (F n ) = {A € C : |A| < 1, A / 1}. 



(10) 



(11) 



4. Finally, since the spectrum of F n is closed, is also contains the point A = 1. Since 
this point is neither in the point spectrum nor in the residual spectrum, it must be 
in the continuous spectrum, that is 



a c (F n ) = {1}. 



(12) 



We proceed with the detailed analysis of each step. 
Step 1. The equation Q can be written as 

= Axi - x n+ i - x n+2 - x n+3 
xi = Ax 2 , 
x 2 = Ax 3 , 

x k = Axfc+i, 



(13) 



Since A = implies x = 0, zero is not an element of a p (F n ). If A ^ 0, by applying ([13 
recursively, we have 



1 



1 



1 



Xk+l = ^ Xk = J2 Xk ~ l = ~\3 Xk - 2 



k > 1. 



Thus 
and 



A* I 



\X 1 = Xl 



y— 

^ \\\ k ' 



(14) 

(15) 
(16) 



If |A| < 1, then ||x||i = oo, so x ^ I 1 . If |A| > 1, then ||x||i = |xi| |A|/(|A| — 1). Inserting 
(I14p into the first equality of fjl3[) gives 



= Axi - x n+ i - x n+2 - x n 

= \ X1 - — X1 -— X1 



+3 



A n+2 



Xi 



Xi 



1 A, 1 1 1 
a -a4 1+ a + a^ + --- 



xi A 



Xl- 



1 1 



1 



A n+1 - A w 
A™- 1 (A - 1) 
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Since x\ 7^ 0, we conclude that a p (F n ) consists of those roots of the polynomial 

Pn+i(A) = X n+1 -X n -l (17) 

for which |A| > 1, as stated in (JTJ)1U 

Since p n+ i(l) = — 1 < and p' n , 1 (A) > for A G R, A > 1, that is, p n +i is strictly 
increasing for A > 1, we conclude that F n has exactly one real eigenvalue larger than one. 
Let us denote this eigenvalue by \ max (F n ). By Ostrovsky's theorem [Pra041 Theorem 
1.1.4, p. 3], A max (F n ) is the unique positive root of p n +i(X) and the absolute values of all 
other roots are strictly smaller. Consequently, all other eigenvalues of F n are in absolute 
value strictly smaller than A max (^n) which, in turn, implies 

ra(F n ) = A max (F n ). (18) 

Figure [U shows cr p (F n ) for various values of n. 
Step 2. The equation ([8]) can be written as 

Vi = Axi - x n+ i - x n+2 ~ %n+3 , 

X2 = j (xi +V2), 

£3 = ^(^2+2/3), (19) 



Xk+i = j (xk + Vk+i), 



By setting 



and using (fl9|) . we have 



00 00 
u= ^2 x k , v = ^2 

k=n+l k=n+l 



111 

u = -x n + -u + -v, yi = Xx 1 -u. 
After rearranging, we have 

1 r 

U = J—^{Xn + V). 

Thus, 

yi = Xxi - j— -{x n + v). (20) 



1 These roots are the eigenvalues and the vectors x defined by (|15p are the corresponding eigenvectors. 
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-1 -0.5 0.5 1 1.5 



-1 -0.5 0.5 1 1.5 



n=25 



n=50 





-1 -0.5 0.5 1 1.5 



-1 -0.5 0.5 1 1.5 



Figure 1: The point spectra cr p (F n ) for various values of n. 



By recursively applying (|19p . we have 



1 



X2 = -r (xi +2/2), 

1 1 1 1 

2^3 = J (X2 + y 3 ) = -2 Xl + ^ y 2 + V 2/3, 



1 . ,11 1 1 1 

£fc+i = j ( x k + Vk+i) = p: xi + y 2 + 2/3 + -3^2 2/4 H ^ y 2/fc+l> ( 21 ) 



A fc -1 ' x h»* ' \k- 



A fc " 



A 



Inserting x n into (|20|) gives 



1 / 1 



yi = Xxi ~ yrj { x^T Xl + a^T y2 + a^2 2/3 + • • • + t Vn + v ) , 



1 

A 
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and solving for x\ gives 
1 

Xl 



X n -\X - l) yi + y 2 + Xy 3 + X'y 4 + • • • + \ n ~ z y n + X 



n-2 , 



\n+l _ \n _ 1 
By inserting this into (|2ip we have 



v n— 1 



{\I-F n )- l y 



\n+l _ A™ - 1 

where the matrix representations of A and B are given b}0 



(A + £)y 



/ (A - l)A n ~ 1 
(A - l)\ n - 2 

(A-l)A 
(A-l) 
(A-l)i 



1 


A 


A 2 


A 3 




1 

A 


1 


A 


A 2 




1 


l 


l 




i 


A n-2 


A n-3 


A n-4 




? 


A n-1 


X"-' 2 


A n-3 




Y 




X n-1 


A n-2 







A n-2 
A n-3 



1 

j_ 

A 
1 

A? 



A n-1 A n-1 
A n— 2 A n —2 



A 
1 

A 



A 

1 

l 

A 



and 



B 



/ 

i 

o X t o o 

Y i i 

U P" A 7 A 

V i i : 







respectively. Obviously, for |A| > 1 we have ||A||i < oo and ||-B||i < oo. Thus, for |A| > 1 
and A not being the root of A n+1 — A™ — 1, the operator XI — F n has a bounded inverse, 
so the resolvent set of F n is given by (|9|). 

Step 3. The point spectrum of the transposed operator consists of all A € M such 
that 

(AJ-F*)x = 0, x t^O, ||x||oo < oo. 

This is equivalent to 

X2 = Axi, 

x 3 = \ x 2 = A 2 Xl, 



x n+l — Ax n — X n Xl, 

x n+2 = Xx n+ i - xi = (X n+1 - 1) Xl, 
x n+3 = Xx n+2 - xi = (X n+2 - X - 1) Xl, 



x k = {X k 



X - l)x 1 , 



Next row of A is obtained by dividing the previous row by A. 
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Therefore, 

/ \k— n— 1 _ i ' 
Xfc = A*" 1 I Xl . 



A-l 

For |A| < 1, A / 1 we have 

kfcl < U + | A 2 _ X | ) i-'-i I- 

which implies ||x||oo < oo. For A = 1 we have 

x 2 = xi, 
x 3 = xi, 



X n +2 = 0, 

= — 2 JCl, 



Xfc = —(k — n — 2) xi, 



so II^Hoo = oo. We conclude that the point spectrum of is given by (ITU]) . This, in turn, 
implies ([TT]) and (fT2]) as described before. 



2.1 Relationship to generalized Fibonacci sequences 

In this section we describe the relationship between operators F n and generalized Fibonacci 
sequences. A generalized Fibonacci sequence {f^} is defined by 

fi H) = h i n) = !,-••, /a = l, ft ] = f^i + ft2 n -v k>n + l. (22) 
For n = 1 this definition yields the classical Fibonacci sequence 

/i = l, / 2 = 1, fk = fk-l + fk-2, k>2. (23) 



7 



By induction we can prove that the £;-th power of the matrix F n from ([5]) for k > n has 
the form 

\ 



An) 
J k—n 


An) 

J k-n+1 


An) 

J k-n+2 


An) 
Jk-1 


An) 
■Ik 


An) 
Jk 


An) 
Jk 


An) 

Jk-n-1 


An) 
J k—n 


An) 

J k-n+1 


An) 
Jk-2 


An) 
Jk-1 


An) 
Jk-1 


An) 
Jk-1 




An) 
J2 


An) 
J3 


An) 
Jn 


An) 
Jn+1 


An) 
Jn+1 


An) 
Jn+1 





1 


1 


1 


1 


1 


1 








1 


1 


1 


1 


1 





1 










1 



1 







1 







We conclude that 



iK fc iii = i + E^ 



i=l 

By applying (|22p to the terms in parentheses we have 
k 



i=l 



2 E ^ = fl U) + • • • + fn^ + (fSl + fnt + • • • + ft\) + f; 



An) 



+ (fi n) + ft' + • • • + f£!n-l) + f?2 n + f^n+1 +- + /, 

f[ n) + f ( 2 n) + --- + fn n) + & + ■■■ + 
, An) , An) An) An) 
"T Jk "T •/ fc-n "T •/ fc-n+1 ' "'"•'fc 



c (n) 



c (n) 



c (ra) 



(n) 



(24) 



EAn) _ An) , An) An) An) , f 

Ji "T Jk-n^ Jk-n+l' " ' ~r Jk ' J i 



e (n) 



(n) 



A >n , An) 



i=l 



Prom this, by applying (|22p again recursively, we obtain 

k 



i=l 



(») 



„(n) „(n) 
■/ k-n ' J k-n+1 



+ --- + f ( k n) +ft ] 



An) , An) An) _ . 

Jk — n+1 ' Jk ^ Jk+1 1 



An) 

'k-n+2 



+ 



An) , f 
Jk+1 ' Jk+2 



in) 



Inserting this into ()24j) gives 



An) 

Jk+n+l 



\\F% 



An) 

Jk+n+l 



(25) 
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and from (1181) it follows that 



k-^oo 



Also, by using standard techniques in analyzing linear recurrence relations with constant 
coefficients, we can prove that for all i,j'H 



lim [F " ]m ' 



lim = X max (F n ) 



m— >oo f( n ) 
J Jra 



For example, by setting n = 1 we have for the Fibonacci sequence (I23p 

lim (A +2 ) 1/fc =r,(F 1 ) = ^, 

fc— >-oo / 

l im A±l = 1 + ^ 

k-^oo ff, 2 

3 Fibonacci-like operators 

Now we would like to consider the family of linear operators r n : I 1 — > I 1 defined by 



P ^ (k- n)x k ,xi,x 2 ,x 3 , ■ ■ • J , n = 1,2,3, 

k=n+l I 



for some real positive The domain of T n is 

Dom r n = jx G Z 1 : 
and its matrix representation in the standard basis is 



(26) 



^ (k - n)x k 

k=n+l 



< OO 



( • • • p 2p 3p 4p 5p 

100000 

010000 

000100 

000010 



(27) 



However, the operator T n is not closed as illustrated by the following example. 



3 The proof is derived using the fact that / ; ' n ' has the form /j = aAj nax (_F n ) + ^" =1 aiA|, where 
Xmax{F n ) and Ai are the roots of the characteristic polynomial (fT7)) . and |A max (_F, I )| > |A<|. 
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Example 1 Let us define the sequence {x^} of vectors in I 1 by 

m+n— 1 



x (m) 

Then 



••• — 

m 



x {m) -> (0 •••)*, 
while 

r n x^) = ( P o ••• o £ o •••) t ^(p o, o -..)*. 

Although the point spectrum of T n is defined and can be computed in a standard 
manner (see later) , the resolvent set of T n is empty, which makes the analysis of T n less 
interesting. Instead, we shall consider the family of operators G n : l\ — > l\ formally defined 
by 

G n — D n V n D n , 

where 

n 

D n = diag(l,--- , l,p,2p,3p,4/v). 
That is, for n E N the operator G n is defined by 

. ( ^ 3 4 5 \ 

(Xi,X 2 ,X 3 , ■■■) ->■ I 2^ ^,£1,^2, • • • ,^n-l,/52;n,2x n+ i,-a; n+ 2,-X n+ 3,-X n+ 4,--- I , 
\fc=n+l ~ ° / 

and its matrix representation in the standard basis is 






•• 


• 


1 


1 


1 


1 


1 


1 


•• 


• 




















1 •• 


• 




















•• 


• p 




















•• 


■ 


2 

















•• 


• 





3 
2 














•• 


• 








■1 

3 











•• 


• 











5 

1 






V 



Let us define the polynomial q n+ \(X) by 

q n+1 (X) = X n+1 -2X n + X n - 1 



(28) 



(29) 



Similarly as in section [2j the spectrum of G n is classified in several steps which are sum- 
marized as follows: 
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1. first, by solving the equation 

(XI -G n )x = 0, x / 0, (30) 
we show that the point spectrum is 

cr P (G n ) = {X G C : e?„ + i(A) = 0, |A| > 1}, n > 2. (31) 

2. second, by solving the equation 

(XI-G n )x = y, x^O, (32) 

we compute the inverse (XI — G n ) _1 and show that the resolvent set consists of all 
A such that |A| > 1 which are not in a p (G n ), 

p(G n ) = {X G C : |A| > 1, A i a p (G n )}, (33) 

3. third, we analyze the transposed operator G l n and show that 

a p (Gi) = {X G C : |A| < 1, A + 1}, (34) 
which, together with ([2]), implies that the residual spectrum of G n is 

a r (G n ) = {A G C : |A| < 1, A + 1}. (35) 

4. Finally, since the spectrum of G n is closed, is also contains the point A = 1. Since 
this point is neither in the point spectrum nor in the residual spectrum, it must be 
in the continuous spectrum, that is 

a c (G n ) = {1}. (36) 

The proofs are similar to the ones from section [21 but more tedious. 
Step 1. The equation (|30p can be written as 

= Axi - x n+ i - x n+2 - x n+3 , (37) 

x k = Xx k+ i, k = !,-■■ ,n-l, 
p x n = X x n +i, 
k — n + 1 



k — n 



x k = Xxk+i, k = n + l,n + 2, 



Since A = implies x = 0, zero is not an element of a p (G n ). If A ^ 0, by applying (f3"T|) 
recursively, we obtain 

x k = ^k—i k = 2,3, ■ ■ ■ ,n, 

k — n 

x k = P ^ fc _x xi, k = n + l,n + 2, ■■■ . (38) 



11 



This, in turn, implies 




l l 
I 5? 



1 JL 



2p 3p 



(k—n) p 



so that 



rc 1 1 1 = xi 




If |A| < 1, then ||x||i = oo, so x ^ i . If |A| > 1, then, by using differentiation of the 
geometric series, we have 



thus, x E I . By inserting (f38j) into the first equality of (i37|) and using differentiation of 
the geometric series we have 



Finally, solving this equation with x\ ^ gives (f3Tj) . 

We shall now prove that (J p {G n ) consists of A max (G n ), a unique simple real eigenvalue 
larger than one and all other eigenvalues have modulus smaller than A max (G n ). This also 
implies 



The proof is based on the ideas from the proof of |Pra04l Theorem 1.1.4, pp. 3-4]. Indeed, 



Since q n+ \{\) = —p < and q' n+ i(\) > for A G R, A > 1, that is, q n+ i(X) is strictly 
increasing for A > 1, we conclude that q n+ i(X) has exactly one real root larger than one 
or, equivalently, that G n has exactly one real eigenvalue larger than one. Let us denote 
this eigenvalue by A max (G n ). Let z / A max (G n ) be some other real or complex eigenvalue 
of G n and let £ = \z\ > 1. Since z is also the root of q n+ i(X), the relation (f40l) implies 






(39) 




z n -\z-l) 2 =p, 



which, in turn, implies 



z 



n-l 



z - 1| 2 = p. 



(42) 



Since £ > 1, this implies 



C--\t-if< P , 
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or 

gn + i(C) = C^ 1 (C-l) 2 -p<0. (43) 

Since q n +%(X) is strictly increasing for A > 1, and q n +i(X max (G n )) = 0, we conclude that 
C — A max (G n ) and that the equality in (|4"3j) holds only if ( = A max (G n ). But, the equality 
in (j43]) and ([32]) imply 

\z-l\=(-l, 

that is, z G R and z = ±( = ±A max (G n ). The choice z = — A max (G n ) is impossible since 
q n+ i{— A max (G n )) ^ 0, and the second choice contradicts the assumption z ^ A max (G n ). 
Therefore, Q < A max (G n ) as desired. 

Remark 1 Although the above analysis is sufficient for our purposes, by inspecting the 
polynomial g n +i(A) and its derivative from (|40p and (jUJ), respectively, we can establish 
further facts about its roots. From (flll) we see that the derivative q' n+1 (X) has exactly 
two real positive simple roots Ai = and A2 = 1 and, if n > 2, also the root Ao = 0. If 
n > 3 then Ao is multiple. Let po = 4A" _1 /(n + l) 2 . The number of real roots of g ra+ i(A) 
in the open interval (0, A max (G n )) is governed by p as follows: if p > po, then there are no 
such roots, if p = po there is exactly one root equal to Ai and if p < po there are exactly 
two roots, one smaller and one larger than Ai. Finally, if n is odd, then q n+ i(X) also has 
a simple negative real root. As we have already proved, A max (G n ) is the root with strictly 
maximal absolute value. 

Remark 2 It is easy to see that the point spectrum of T n from (|26|) and (|27p is equal to 
the point spectrum of G n . 

Step 2. The equation (1321) can be written as 

yi = Axi - x n+1 - x n+2 - x n+3 , (44) 

Xk+i = -r (xk + Vk+i), k = 1, • • • ,n - 1, 

x n +i = -r (px n + y n +l), 

1 (k-n + l \ 
Xk+i = T ( k _ n x k + Vk+i J, k = n + 1, n + 2, ■ ■ ■ . 

By recursively applying the above three equalities, we obtain 

(k-n)p (k-n)p (k - n) p (k - n) p 

X *= X k-1 X l+ X k-1 X k-2 S/3 + ---+ Afc _ TO+1 Vn+ (45) 

k — n 1 k — n 1 k — n 1 

H ^ -JkZ^ Vn+l H ^ X k-n-l yn + 2 ^ 3 \k-n-2 y ™+ 3 ^ 

k — n 1 1 
+ k - n - 1 A 2 Vk ~ X + A Vk ' k = n + l,n + 2,--- . 

For |A| > 1, by inserting (|45[) into (|44p . rearranging, and using differentiation of the 
geometric series, we have 

a2 ( P P P P , 1 



f 1 = A 001 ~ (a -1)2 ^ + >7 y2 + A^T ys + ' ' ' + A 2 yn + A Vn+1 

- «n+2 Un+2 ~ «n+3 yn+3 ~ «n+4 2/n+4 , (46) 
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where 



1 fc + 1 1 k + 2 1 k + 3 1 
an+k ~ A + ~k~ A 2 + ~k~ A3 + ~T~ A* + 



= 2,3,4, 



and 



\a n +k\ < 



1 / ^k + i 1 \ 1 / IN 



|A| 



Solving for (j36]) for xi gives 
1 



A n - 2 (A-l) 2 yi + py 2 + pAy 3 + pA 2 y 4 + • • • + p\ n ~ l Vn + X n ~ L y n . 



(|A|-1) S 



n— 2 „, i \n— 1 



9n+l(A) 

+ 0>n+2 Vn+2 + OL n+ z Un+3 + «n+4 Vn+4 + 



where 



and 



X' 



\a n +k\ < |A| 
By inserting (j3HJ) into ([34]) and ([15]) . we have 

x = (A/-G n )- 1 y = 



2 (A - l) 2 a n+ k, 

n-1 (|A| + 1) 2 



(|A|-1) 



2 ' 



1 



2«+i(A) 

where the matrix representations of A and B are given by 

pA 



(A + B)y 



( (A 
(A 
(A 
(A 



1) 2 A"" 2 

1)2A n-3 
l)2 A n-4 
!)2 A n-5 



(A-1) 2 A 
(A -I) 2 
(A"1) 2 A 
(A-l) 2 ^ 
(A-l) 2 ^ 
(A-l) 2 f£ 
(A-l) 2 ^ 



P 

£ 

A 
P 

A? 
P 

F 



p 

P 
P 

eL 

A n 

2P 2 
3P 2 
4p 2 



P 

P 
A 
P 



A™^ 
P 

A™^ 
P 

P 2 
3P 2 

A n+1 
4p 2 



PA 2 
pA 



p 

A 



A n-5 
P 

A™~^ 
P 

P 2 

2p 2 

3p! 
A" 
4P 2 



+ 3 A n + 2 A™ +1 
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respectively. For |A| > 1, by using differentiation of the geometric series, ([49)) . and the 
argument used in (|4"T|) . we have < oo and < oo. Thus, for |A| > 1 and 

q n +i(X) 0, the operator XI — G n has a bounded inverse and its resolvent set is given by 

Step 3. The point spectrum of the transposed operator G l n consists of all A £ M, 
|A| < 1, such that 

{XI — G t n )x = 0, x / 0, \\x\loo < oo. 

This is equivalent to 



x k = X k 1 xi, 



,n, 



P 

*„ +fc = i ( - A fc - 2 - 2 A fe - 3 - 3 A fc ~ 4 

If A ^ 1, we have 



(k-2)X-(k-l))x 1> k>2 



1 

•^n+k ~j~ 



^n+k- 
P 



(X k - 2 + 2A 



fc-3 



+ ••• + 



2) A +(*-!)) 



(A -I) 2 
(A -I) 2 J 



it 



P 



Therefore 



(A -1)2 
1/1 l + k + k-l 



lXn+kl ~k{p + |A-1| 2 



x\ < 



xi, fc > 2. 



+ 



xi, 



%n+k 



1 



which implies ||x||oo < oo. For A = 1 we have 

-(- 

k \p 

1/1 (jfe-l)jfc 
k\~p 2 



Xi, 



2p |A-lf 



{k-2)-{k-l))x x 
k > 2, 



k > 2, 
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so 1 1 x I |oo = oo. We conclude that the point spectrum of G l n is given by (f34"|) . This, in turn, 
implies (j35l) and (j36j) as described before. 



Example 2 The lesion forming plant pathogen potato late blight (phytophthora infes- 
tans) grows radially on a leaf with a constant daily rate. The latency period for a lesion 
to become infectious is five days, and the sporulating area is infectious for one day. In 
|PSV05j the epidemic spread of such pathogen is modeled with the infinite dimensional 
Leslie matrix of the form of T$ as defined in (j27f) . Further, the upper bound for the 
speed of invasion in computed via minimization of the largest eigenvalue A max (r5(s)). 
From Remark [2] it follows that this eigenvalue is the largest unique positive root of qe(X) 
from (I29p . Here the parameter p has the form p(s) = const x M{s) where M(s) is some 
moment-generating function (for example, M(s) = exp(cr 2 s 2 /2) for the Gaussian kernel 
or M(s) = 1/(1 — a 2 s 2 ) for the Laplace kernel). Here T^(s) appears naturally due to the 
fact that the considered pathogen has a latency period of five days. It is interesting that 
A m ax(r5(s)) can be computed analytically: 

. (r , _ 1 2V3 (2 + 27Vp+^10Vp + 729p) 1/3 

3 3 (2 + 27^/Ti + VlOS^p + 729p) 1/3 3 • 2Vi 

The speed of invasion is bounded by 

v* = min - In |A max (r 5 (s))] , 

0<s<s S 

where s is the maximum s for which M (s) is defined. For details about a rather complex 
derivation of this model we refer the reader to [PSV05] . 
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